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What is claimed is : 

1. A method for reducing noise in a digital image 
comprising pixels, the method comprising the steps of: 

computing global statistics from the image; 
for each of a plurality of image pixels: 

computing local statistics for the pixel; 
using the local and global statistics to 
configure a local filter for filtering that pixel; and 

filtering the pixel using the local filter 
to reduce image noise. 

2. The method of claim 1, wherein the step of 
computing the global statistics further comprises the steps 
of estimating the global noise standard deviation cr to 
generate the global statistics. 

3. The method of claim 1, wherein the step of 
computing the local statistics for each pixel further 
includes the steps of: 

selecting a window containing said pixel and a 
plurality of neighboring pixels; 

computing the 2-D local variance of said pixel 
based on information related to the pixels in the window; 
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computing the 1-D local variances along multiple 
directions through said pixel within the window; and 

detecting the local edge direction by selecting 
one of the direction with the smallest 1-D local variance. 

5 

4. The method of claim 1, wherein the step of 
computing the local statistics for each pixel further 
includes the steps of: 

selecting a window containing said pixel and a 
10 plurality of neighboring pixels; 

computing the 2-D local variance a] of said pixel 
based on information related to the pixels in the window; 

computing the 1-D local variances a], a\ , a\ , and 

a\ along the horizontal (L x ) , vertical (L 2 ) , diagonal from 
15 upper left to lower right (L 3 ) , and diagonal from upper 
right to lower left (L 4 ) directions through said pixel, 
respectively, within the window; and 

detecting the local edge direction by selecting 
one of the directions with the smallest 1-D local variance. 

20 

5. The method of claim 1, wherein the step of 
configuring the local filter for each pixel using the local 
and global statistics further includes the steps of: 
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selecting the detected local edge direction L as 
the direction of the local filter; 

for the detected local edge direction L computing 
the 1-D filter strength as a function of the square root of 
5 the local variance and the global noise standard deviation; 

computing the 2-D filter strength as a function 
of the local variance and the global noise standard 
deviation; and 

configuring the local filter for the detected 
10 local edge direction L based on the 1-D and 2-D filter 
strengths . 



6. The method of claim 5, wherein the step of 
configuring the local filter for each pixel using the local 
15 and global statistics further includes the steps of: 

selecting the detected local edge direction L k (k 
= 1, 2, 3, or 4) as the direction of the local filter; 

for the detected local edge direction L k computing 
the 1-D filter strength a k = min(2<7,max(3<7 - cr^ ,0)) /(2<r) ; 
20 computing the 2-D filter strength 

a 0 = min(2cr, max(3<7 - a 0 ,0)) /(2a) ; and 
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configuring the local filter f k for the detected 
local edge direction L k according to the following 
conditions : 

(i) if the detected direction is Li, then f x is 
configured as a 2-D local filter for horizontal direction, 
wherein: 



a n 



a 0 +3a,(l-a 0 ) a 0 +3(3-2a,)(l-a 0 ) a 0 + 3a,(l-a 0 ) 



(ii) if the detected direction is L2, then f 2 is 
configured as a 2-D local filter for vertical direction, 
wherein : 



/ 2 = 



1 



(iii) 



a 0 a 0 +3a 2 (\-a 0 ) a 0 ' 
a 0 a 0 +3(3-2or 2 )(l-a 0 ) a Q 
a 0 a 0 +3a 2 (l-a 0 ) 



if the detected direction is L 3 , then / 3 



is configured as a 2-D local filter for the diagonal 
direction from upper left to lower right, wherein: 



9 



a 0 +3a 3 (l-« 0 ) 



a* 



a 0 + 3(3-2a 3 )(l-a 0 ) 



a. 



a 0 + 3« 3 (l-a 0 ) 



and 
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(iv) if the detected direction is L 4 , then f 4 is 
configured as a 2-D local filter for the diagonal direction 
from upper right to lower left, wherein: 



7. The method of claim 2, wherein the steps of 
estimating the global noise standard deviation a further 
includes the steps of: 

dividing the input image into overlapping or non- 
10 overlapping blocks; 

computing the mean and the standard deviation for 
each block; 

finding the smallest standard deviation d 0 and its 
corresponding meanm 0 ; 
15 detecting block saturation due to noise; 

compensating for the smallest standard deviation^ 

to generate a compensated smallest standard deviation d 0 ; 

selecting the block standard deviations d n that 
are within a range of the compensated smallest standard 
20 deviation d 0 ; and 




a 0 

a 0 +3(3-2tf 4 )(l-a 0 ) 



a 0 + 3a 4 (l-a 0 ) 



a Q +3a 4 (l-a 0 ) 



5 
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averaging the selected block standard deviations 
d n to generate an estimate of the global noise standard 
deviations . 

5 8. The method of claim 7, wherein the block size 

is7x7 or 5x9pixels. 

9. The method of claim 7, wherein the steps of 
detecting saturation and compensating the smallest standard 
10 deviation further include the steps of determining the 
following: 

defining is an upper pixel value limit UL, a 
lower pixel value limit LL, and mid value M between UL and 
LL, 

15 wherein if the mean m 0 is less than the mid range 

M, and the smallest standard deviation is greater than the 
difference between the mean w 0 and the lower limit LL, then 
saturation has occurred at the lower limit LL, and the 
smallest standard deviation d 0 is compensated by adding 

20 thereto a compensation term that is a function of the 

smallest standard deviation d 0 and said difference between 
the mean m 0 and the lower limit LL, to generate the 
compensated smallest standard deviation d 0 ; 

28 



SAM2.0033 

else if the mid range M is less than the mean 
m 0 , and the smallest standard deviation d 0 is greater than 
the difference between the upper limit UL and the mean m Q , 
then saturation has occurred at the upper limit UL, and the 
smallest standard deviation d Q is compensated by adding 
thereto a compensation term that is a function of the 
smallest standard deviation d Q and said difference between 
the upper limit UL and the mean m Q , to generate the 

compensated smallest standard deviation d 0 ; 

otherwise, no saturation has occurred, wherein 

d 0 =d 0 . 

10. The method of claim 7, wherein the steps of 
detecting saturation and compensating the smallest standard 
deviation further include the steps of determining the 
following: 

where UL is an upper pixel value limit, LL is a 
lower pixel value limit, and UL < M < LL, if the mean 
m 0 <M and the smallest standard deviation d 0 >m Q -LL , then 
saturation has occurred at the lower limit LL, wherein d 0 is 

compensated as d 0 =d 0 +K-(d Q -(m 0 -LL)) , such that K is a 
compensation factor; 
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else if the mean m Q >M and the smallest standard 
deviation d 0 >UL-m 0 , then saturation has occurred at the 
upper limit UL, wherein d 0 is compensated as 

d Q =d 0 +K-(d Q -(UL-m 0 )); 

otherwise, no saturation has occurred, wherein 

d 0 =d Q . 

11. The method of claim 10, wherein LL = 0, UL = 255, 
and M = 128. 

12. The method of claim 7, wherein the step of 
selecting the block standard deviations further includes 
the steps of selecting the block standard deviation d n for 

averaging if \d n -d 0 \< max(t/ 0 ,l) . 

13. A noise reduction system for reducing noise in a 
digital image comprising pixels, the system comprising: 

a global statistics module that computes global 
statistics from the image; 

a local statistics module that computes local 
statistics for each of a plurality of image pixels; 
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a filter configuration module that uses the local 
and global statistics for a pixel to configure a local 
filter for filtering that pixel; and 

a filter that as configured filters the pixel to 
reduce image noise. 

14. The system of claim 13, wherein the global 
statistics module estimates the global noise standard 
deviation a to generate the global statistics. 

15. The system of claim 13, wherein the local 
statistics module computes the local statistics for each 
pixel by: 

selecting a window containing said pixel and a 
plurality of neighboring pixels; 

computing the 2-D local variance of said pixel 
based on information related to the pixels in the window; 

computing the 1-D local variances along multiple 
directions through said pixel within the window; and 

detecting the local edge direction by selecting 
one of the direction with the smallest 1-D local variance. 
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16. The system of claim 13, wherein the local 
statistics module computes the local statistics for each 
pixel by: 

selecting a window containing said pixel and a 
5 plurality of neighboring pixels; 

computing the 2-D local variance a\ of said pixel 
based on information related to the pixels in the window; 

computing the 1-D local variances a] , a\ , a], and 

a] along the horizontal (Li) , vertical (L 2 ) , diagonal from 
10 upper left to lower right (L 3 ) , and diagonal from upper 
right to lower left (L 4 ) directions through said pixel, 
respectively, within the window; and 

detecting the local edge direction by selecting 
one of the directions with the smallest 1-D local variance. 

15 

17. The system of claim 13, wherein the filter 
configuration module configures the local filter for each 
pixel using the local and global statistics by: 

selecting the detected local edge direction L as 
20 the direction of the local filter; 

for the detected local edge direction L computing 
the 1-D filter strength as a function of the sguare root of 
the local variance and the global noise standard deviation; 
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computing the 2-D filter strength as a function 
of the local variance and the global noise standard 
deviation; and 

configuring the local filter for the detected 
local edge direction L based on the 1-D and 2-D filter 
strengths . 

18. The system of claim 17, wherein the filter 
configuration module configured the local filter for each 
pixel using the local and global statistics by: 

selecting the detected local edge direction L k (k 
= 1, 2, 3, or 4) as the direction of the local filter; 

for the detected local edge direction L k computing 
the 1-D filter strength a k = min(2cr, max(3cr - <j k ,0)) /(2a) ; 

computing the 2-D filter strength 
a 0 = min(2cr, max(3cr - a 0 ,0)) /(2a) ; and 

configuring the local filter f k for the detected 
local edge direction L k according to the following 
conditions : 

(i) if the detected direction is Li, then /, is 
configured as a 2-D local filter for horizontal direction, 
wherein: 
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h 9 



a 0 +3a,(l-a 0 ) a 0 + 3(3 - 2ar,)(l - « 0 ) a 0 +3a,(l-a 0 ) 



(ii) if the detected direction is L 2 , then f 2 is 
configured as a 2-D local filter for vertical direction, 
wherein : 



»5 



(iii) 



a 0 a 0 +3or 2 (l-a 0 ) a 0 
a 0 a 0 +3(3-2a 2 )(l-or 0 ) a 0 
a 0 a 0 +3a 2 (\-a 0 ) a 0 



if the detected direction is L 3 , then f } 



is configured as a 2-D local filter for the diagonal 
direction from upper left to lower right, wherein: 



h 9 



a 0 +3« 3 (l-a 0 ) 



a 0 +3(3-2a 3 )(l-a 0 ) 



a 0 +3a 3 (\-a 0 ) 



10 and 



(iv) if the detected direction is L 4 , then / 4 is 



configured as a 2-D local filter for the diagonal direction 
from upper right to lower left, wherein: 



/ 4 = 



1 



a 0 +3or 4 (l-a 0 ) 



a 0 +3(3-2a 4 )(l-a 0 ) 



_a 0 +3a 4 (\-a 0 ) 



15 
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19. The system of claim 14, wherein the input image 
comprises overlapping or non-overlapping blocks, and 
wherein the global statistics module further comprises: 
a mean and standard deviation module that 
5 computes the mean and the standard deviation for each block 

a minimum finder module that finds the smallest 
standard deviation d 0 and its corresponding meanw 0 ; 

a saturation detector that detects block 
10 saturation due to noise; 

a saturation compensator that compensates for the 
smallest standard deviation^ to generate a compensated 

smallest standard deviation d 0 ; and 

a selective averaging module that selects the 
15 block standard deviations d n that are within a range of the 

compensated smallest standard deviation d 0 , and averages 
the selected block standard deviations d n to generate an 
estimate of the global noise standard deviations" . 

20 20. The system of claim 19, wherein the block size 

is7x7 or 5x9pixels. 

21. The system of claim 19, wherein: 
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an upper pixel value limit is denoted UL, a lower 
pixel value limit is denoted LL, and a mid value M is 
between UL and LL, 

wherein the saturation detector determines if the 
5 mean m Q is less than the mid range M, and the smallest 
standard deviation is greater than the difference between 
the mean m 0 and the lower limit LL, indicating that 
saturation has occurred at the lower limit LL, and if so, 
the saturation compensator compensates the smallest 
10 standard deviation d Q is by adding thereto a compensation 
term that is a function of the smallest standard deviation 
d 0 and said difference between the mean m 0 and the lower 
limit LL, to generate the compensated smallest standard 
deviation d 0 ; 

15 else if the saturation detector determines that 

the mid range M is less than the mean m 0 , and the smallest 
standard deviation d 0 is greater than the difference between 
the upper limit UL and the mean m 0 , indicating saturation 
has occurred at the upper limit UL, the saturation 

20 compensator compensates the smallest standard deviation d 0 
by adding thereto a compensation term that is a function of 
the smallest standard deviation d 0 and said difference 
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between the upper limit UL and the mean m 0 , to generate the 

compensated smallest standard deviation d 0 ; 

otherwise, no saturation has occurred, wherein 

d 0 =d 0 . 

22. The system of claim 21, wherein: 

if the mean m 0 < M and the smallest standard 
deviation d 0 >m Q -LL, indicating saturation has occurred at 
the lower limit LL, then d 0 is compensated as 
d 0 =d 0 +K-(d 0 -(m 0 -LL)) , such that ^ is a compensation factor; 

else if the mean m 0 >M and the smallest standard 
deviation d 0 >UL-m 0 , indicating saturation has occurred at 
the upper limit UL, then d 0 is compensated as 

d 0 =d 0 + K-(d 0 -(UL-rn Q )); 

otherwise, no saturation has occurred, wherein 

d Q =d 0 . 

23. The system of claim 22, wherein LL =0, UL = 255, 
and M = 128 . 
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24. The system of claim 19, wherein block standard 
deviations d n are selected for averaging if | d n - d 0 |< max(^ 0 
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